Local electronic structure of Fe 2+ impurities in MgO thin films: Temperature 
dependent soft x-ray absorption spectroscopy study 
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We report on the local electronic structure of Fe impurities in MgO thin films. Using soft x-ray 
absorption spectroscopy (XAS) we verified that the Fe impurities are all in the 2+ valence state. 
The fine details in the line shape of the Fe Li,z edges provide direct evidence for the presence of 
a dynamical Jahn- Teller distortion. We are able to determine the magnitude of the effective Da.h 
crystal field energies. We also observed a strong temperature dependence in the spectra which 
we can attribute to the thermal population of low-lying excited states that are present due to the 
spin-orbit coupling in the Fe 3d. Using this Fe 2+ impurity system as an example, we show that an 
accurate measurement of the orbital moment in Fe3C>4 will provide a direct estimate for the effective 
local low-symmetry crystal fields on the Fe 2+ sites, important for the theoretical modeling of the 
formation of orbital ordering. 



PACS numbers: 71.70.Ch, 71.70.Ej, 75.10.Dg, 78. 70. Dm 



Magnetite is one of the most controversially discussed 
systems in solid state physics. 1 It shows a first order 
anomaly in the temperature dependence of the electrical 
conductivity at 120 K, i.e., the famous Verwey transi- 
tion^ which is accompanied by a structural phase tran- 
sition from the cubic inverse spinel to a distorted struc- 
ture. It is only very recently that one realizes that this 
transition may involve not only charge ordering of Fe 2+ 
and Fe 3+ ions but also t2 g orbital ordering at the Fe 2+ 
sites^— Important in this regard are the recent results 
from band theory studies^ in which the charge and or- 
bital occupations were calculated based on the available 
crystal structure data£i 

It is highly desired to determine experimentally the 
electronic structure of Fe3C>4 and especially the local en- 
ergetics of the Fe 2+ sites in order to test the conditions 
under which the ti g orbital polarization and ordering can 
occur. Unfortunately, a direct approach to this system 
is difficult since the simultaneous presence of Fe 2+ and 
Fe 3+ valences as well as octahedral and tetrahedral sites 
makes standard electron spectroscopic methods to yield 
rather broad spectral line shapes, i.e., too featureless for 
a precise analysis concerning the details about the effec- 
tive crystal fields with a symmetry lower than Oh& 

Here we report on our study of the electronic struc- 
ture of Fe impurities in MgO thin films. Having the lo- 
cal quasi-octahedral (Oh) symmetry and similar metal- 
oxygen bond lengths, the impurity system could serve as 
a valuable reference for the more complex Fe 2+ contain- 
ing magnetite. Using soft x-ray absorption spectroscopy 
(XAS) and an analysis based on full multiplet cluster cal- 
culations we found that the Fe impurities are all in the 
2+ charge state and that a dynamical Jahn- Teller dis- 
tortion is clearly present. The spectra showed a strong 



temperature dependence which can be traced back to the 
existence of low-lying excited states due to the presence 
of the spin-orbit interaction. We were able to make es- 
timates concerning the magnitude and temperature de- 
pendence of the orbital and spin contributions to the lo- 
cal magnetic moments. We infer that these local effects 
need to be included when interpreting the temperature 
dependence of the orbital and spin moments in magnetite 
across the Verwey transition. 

Fe^Mgi.^O samples were prepared as polycrystalline 
thin films in an ultra-high vacuum molecular beam 
epitaxy (MBE) system with a base pressure of 5 x 
1CU 10 mbar. High purity Mg and Fe metal were co- 
evaporated from alumina crucibles onto clean Cu sub- 
strates. Molecular oxygen was simultaneously supplied 
through a leak valve. The oxygen partial pressure 
was kept at about 1 x 1CU 7 mbar and monitored by a 
quadrupole mass spectrometer during growth. The Mg 
effusion cell temperature was kept at 315 °C correspond- 
ing to a Mg deposition rate of 2.6 A/min as verified using 
a quartz crystal thickness monitor. During growth the 
substrate temperature was kept at 250 °C which led to 
a distillation process that allows the growth of stoichio- 
metric MgO. This resulted in a MgO deposition rate of 
about 4 A/min. The thickness of the Fe x Mg 1 _ 2 ,0 films 
was about 120 A. The use of thin films on metallic sub- 
strates was necessary to avoid charging problems which 
otherwise readily occur in electron spectroscopic experi- 
ments on bulk insulators like MgO. Clean Cu substrates 
were prepared by growing roughly 1000 A thick Cu films 
in-situ on top of atomically-flat, epi-polished and oxygen- 
annealed MgO substrates. 

The XAS measurements were performed at the 11 A 
Dragon beamline of the National Synchrotron Radiation 
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FIG. 1. Fe L 2 , 3 XAS spectra of Fe^Mgi-^O films for 
x — 0.04, 0.02, 0.01. Included is also the spectrum of bulk 
Fei-aO (top curve, reproduced from Ref. I 111 ). The under- 
lying dashed lines represent the XAS signal of a pure MgO 
film in the photon energy region of the Fe 1/2,3, scaled to fit 
the pre-edge background of the respective Fe^Mgi-zO spec- 
trum. All spectra were taken at room temperature. 



Research Center (NSRRC) in Taiwan^ The photon en- 
ergy resolution at the Fe L 2 ^ edges {hv « 700-730 eV) 
was set at about 0.35 eV. The spectra were recorded us- 
ing the total electron yield (TEY) method in the normal 
light incidence. The base pressure of the XAS chamber 
was 2 x 10~ 10 mbar. The MBE system was directly 
connected to this XAS chamber so that the freshly pre- 
pared samples could be transferred and measured all in- 
vacuo, thereby assuring the cleanliness and reliability of 
the spectra presented here. 

In Fig. [T]we show the experimental Fe £2,3 XAS spec- 
tra of Fe^Mgi.^O films for x = 0.04, 0.02, and 0.01. We 
have also included the spectra of a bulk Fei_aO crystal 
(reproduced from Ref. 111! ) and of a pure MgO film (un- 
derlying dashed lines) as references. The MgO spectra 
have been scaled to fit the pre-edge background of the re- 
spective Fe^Mgi-^O spectrum. The pure MgO film was 
prepared under the same conditions as the Fe^Mgi-^O 
samples, i.e., it has x = 0.00. The Fe £2,3 spectra are 
dominated by the Fe 2p core-hole spin-orbit coupling, 
which splits the spectrum roughly into two parts, namely 
the L 3 {hv k, 708 eV) and L 2 {hv w 721 eV) white line 
regions. The line shapes of the spectra depend strongly 
on the multiplet structure given by the atomic-like Fe 2p- 
3d and 3d-3d Coulomb and exchange interactions, as well 
as by the surrounding solid. 

In going from bulk Fei^O to the films with decreasing 
Fe concentrations, we can clearly observe that the spec- 
tral features become sharper. This can be taken as an 
indication for the presence of inter-Fe interactions in the 
more concentrated systems. Here we would like to note 
that the presence of Fe 3+ species may also contribute 
to the broad spectral features of bulk Fei-^O, a mate- 



rial known to have inherent defects! 12 ' 13 For Fe concen- 
trations lower than 2% the spectra do not significantly 
change anymore - apparently here we already arrived at 
the impurity limit. We notice that for the lowest Fe con- 
centrations the pre-edge XAS background is increasing. 
This can be attributed to the contribution of the MgO 
to the XAS signal in the Fe £2,3 region as shown by the 
spectrum of the pure MgO film. 

We now focus on the temperature dependence of the 
spectra. Fig. [5] shows a close-up of the experimental 
Fe L 3 and L 2 XAS spectra of Feo.02Mgo.gsO for various 
temperatures ranging from 77 up to 500 K. For clarity, 
we here subtracted the XAS background coming from 
the pure MgO film in the Fe L 2 ,3 region. Clear and sys- 
tematic changes with temperature can be observed in the 
spectra. This can be taken as a direct indication for the 
presence of local low-lying excited states. 

To interpret and understand the spectra and their 
temperature dependence, we have performed simulations 
of the atomic-like 2p 6 3d" -> 2p 5 3d™ +1 (n = 6 for 
Fe 2+ ) transitions using the well-proven configuration- 
interaction cluster model.— ~— Within this method we 
have treated the Fe impurity site as an FeOg cluster 
which includes the full atomic multiplet theory and the 
local effects of the solid. It accounts for the intra-atomic 
3d-3d and 2p-3d Coulomb interactions, the atomic 2p and 
3d spin-orbit couplings, the local crystal field, and the 
O 2p-Fe 3d hybridization. This hybridization is taken 
into account by adding the 3d n+1 L and 3d n+2 L 2 etc. 
states to the starting 3d" configuration, where L denotes 
a hole in the O p ligands. Parameters for the multi- 
pole part of the Coulomb interactions were given by the 
Hartree-Fock values f 14 while the monopole parts {Udd, 
U p d) were estimated from photoemission experiments on 
FeOJi The one-electron parameters such as the O 2p- 
Fe 3d charge-transfer energies and integrals as well as the 
crystal field values were tuned to find the best match to 
the experimental spectra. The simulations were carried 
out using the program XTLS 8.3^^ 

Starting with the simple crystal field scheme, an Fe ion 
in the Oh coordination will have its 3d states split into 
the lower lying t 2g and higher e g levels, with the split- 
ting given by lODq of order 1 eV (see Fig. [3]). For an 
Fe 2+ ion, 5 electrons will occupy all the available spin-up 
states. The remaining electron occupies one of the three 
spin-down t 2g orbitals {t 2g e g ), giving the high spin 5 = 2 
Hund's rule ground state. In a multiplet scheme, neglect- 
ing the spin-orbit interaction, the ground state is termed 
the 5 T 2g state, well separated by about 1 eV or more from 
the higher lying b E g and other lower spin configurations. 
Assigning a pseudo orbital momentum of L = 1 to the 
open t 2g shell ) 19 ' 20 the spin-orbit interaction will couple 
it to the S = 2 spin, resulting in three different states 
with J = 1, 2, and 3. The presence of the spin-orbit cou- 
pling in the 3d shell thus splits the 15-fold 5 T 2g state into 
the J states with degeneracies of 3, 5, and 7, respectively 
(see Fig. Using typical parameters for FeO ) 14 ' 18 we 
find an energy separation between them of about 33 meV 
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FIG. 2. (Color online) Temperature dependence of the experimental Fe L3 and L2 XAS spectra of Feo.02Mgo.9gO after 
subtraction of the pure MgO film background. 
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FIG. 3. (Color online) Energy level diagrams for a Fe 2+ clus- 
ter in the Oh coordination in a crystal field (left) and a full 
multiplet (right) scheme. HS denotes the Hund's rule high 
spin configuration 
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FIG. 4. Theoretical Fe Z/2,3 XAS spectra starting from the 
three lowest multiplet states of the FeOe cluster in Oh sym- 
metry, namely the J = 1, 2, and 3 manifolds of the T^g. 



(~ 383 K) and 51 meV (~ 597 K), respectively. 

Important for the understanding of the line shape of 
the Fe £2,3 XAS spectra and their temperature depen- 
dence is that initial states with different quantum num- 
bers could produce quite different XAS spectra, since the 
dipole selection rules, e.g., A J = 0, ±1, will dictate which 
of the possible final states can be reached in the photo- 
absorption process. This is shown in Fig. 21 Indeed, 
each of the three different J = 1, 2, and 3 states has its 
own characteristic XAS spectrum. It is then also quite 
natural to expect a strong temperature dependence for 
the XAS spectrum of a Fe 2+ system if an increase in tem- 
perature causes a thermal population of the J = 2 and 
3 excited states at the expense of a depopulation of the 
J = I ground state. 

Yet, a detailed comparison between the experimental 
spectra and the simulations for the Oh case reveals im- 
portant quantitative discrepancies. A closer look is pro- 



vided in panel (a) of Fig. [5] One can clearly observe 
that neither the simulated K spectrum, i.e., from the 
pure J = 1 ground state, nor the simulated 77 K spec- 
trum, i.e., containing some amount of the J = 2 excited 
state, can reproduce the experimentally obtained 77 K 
spectrum. In particular, feature A is a single peak in 
the experiment while the Oh simulation produces two 
peaks, and feature B of the experiment has considerably 
more weight than can be generated by the simulation. 
All these strongly suggest that the symmetry must be 
lower than Oh, in line with earlier studies using optical 
and Mossbauer spectroscopies^ - — 

We have investigated two further scenarios: the D^d 
(trigonal) and D±h (tetragonal) cases. We find that the 
D 3 d scenario does not provide a better fit, e.g., feature A 
still has a clear two peak structure and peak B has also 
not enough weight, as depicted in panel (c) of Fig. [5] 
On the other hand, we have been able to obtain a very 



4 





t\ 

A 


, . . . . | . . . . | . . . . | . . 

i Simulation K 

\ Simulation 77 K r\ 












t\ 
A 


\ Simulation K 

Simulation 77 K /v 




V^f^o,* ® 77 ^^^^ ; 




, i , . , . i , . , , i , , , , i , , 




A 

A 


Simulation K 

\ Simulation 77 K f\ 

^ 




V^o^o,, ® 77 ^^^^ ; 









705 710 715 720 725 

Photon energy (eV) 



FIG. 5. (Color online) Comparison between the experimen- 
tal 77 K Fe Z/2,3 XAS spectrum of Feo.02Mgo.9gO with the 
simulated K and 77 K spectra of the Fe06 cluster in the 
(a) Oh, (b) D4h, and (c) D^d symmetry. Parameters for the 
simulations are explained in the text. 

good fit using the scenario as shown in the middle 
panel (b) of Fig. [SJ the simulation gives more weight 
for peak B in better agreement with the experiment, and 
the simulated peak A is a more singly peak now as it is 
in the experiment. The overall line shape is thus well 
reproduced. We note that this also provides evidence 
that the Feo.02Mgo.9sO film contains only Fe 2+ ions since 
the simulation has been done for an FeOe cluster having 
the 3d™, 3d n+1 L, and 3d n+2 L 2 configurations with n = 
6. We thus find no indication for the presence of Fe 3+ 
or Fe 1+ ions in our MBE-grown thin film samples, in 
contrast to other earlier studies . 23 ' 24 ' 26 " — 

Continuing now with the D^h scenario, we also have 
simulated the temperature dependence of the Fe L 3 and 
L2 XAS spectra. Fig. [5] shows the results. One can ob- 
serve that the experimentally obtained temperature de- 
pendence (see Fig. [3]) is quantitatively very well repro- 
duced. The good agreement between simulation and ex- 
periment at all the temperatures measured can be taken 
as a strong indication that the D±h scenario describes ac- 
curately the local symmetry of the Fe ion in MgO and 



that the model parameters chosen are realistic giving also 
an appropriate energy separation between the ground 
state and the excited states. In the following we will 
discuss in more detail the total energy level diagram of 
the Fe 2+ cluster in D^h symmetry. 

As already mentioned above, in Oh symmetry the Fe 
3d 6 5 T2 g ground state is split by the spin-orbit interaction 
into the J = 1, 2, and 3 states. A closer look reveals that 
there are also smaller splittings within the J = 2 and 3 
manifolds. The ground state with J = I is also labeled as 
T^g, while the higher lying first excited states with J = 2 
are given the terms and T^. We now switch on the 
Dih crystal field by introducing the parameters Ds and 
Dtr^ as well as differences in the O 2p-Fe 3d hopping 
integrals along the c-axis vs. the a-axis of the FeOe clus- 
ter. The Dih effective crystal field parameter can then be 
most conveniently described as the effective energy split- 
ting Ai2 9 between the xy and the yz/zx states. This 
splitting can be determined from a total energy calcula- 
tion for which the spin-orbit interaction is set to zero. 
Here a negative Ai2 9 means that the xy is the lowest 
state (compressed octahedron). The resulting total en- 
ergy diagram (including spin-orbit interaction) vs. A£2 g 
is plotted in Fig. [7J 

In going from Oh to D^h with increasing A£2 9 , we find 
that the splitting between the T§ g ground state and part 
of the excited states T^g and becomes reduced. It 
decreases from roughly 30 meV for At2 g = to approx- 
imately 10 meV for At2 g = —76 meV, the value with 
which wc find the best simulations for our XAS data. 
See also the inset of Fig. [7] This 10 meV value agrees 
very well with earlier optical and spin relaxation mea- 
surements which have inferred the existence of a state at 
about 100-115 cm -1 1 26 ' 27 In our simulations we need to 
have the splitting reduced from its large cubic value of 
30 meV to this particular 10 meV number in order (1) to 
have sufficient admixing of the and into the pri- 
marily Tsg-like ground state so that feature A becomes 
more like a single peak and feature B gains substantial 
spectral weight as in the experiment (see Fig. [5]), and 
(2) to obtain sufficient thermal population of the excited 
states with increasing temperature in the 77-500 K range 
so that the experimentally observed strong temperature 
dependence is reproduced (see Figs. [2] and [6]). We would 
like to note that the ground state of this 3d 6 system in 
the D^h symmetry is a singlet, as can be seen in the inset 

of Fig. m 

As indicated above, we found At2 g = —76 meV to be 
the optimal number for the effective energy splitting in 
order to achieve the best simulations for the experimental 
data. It is important to realize that this splitting lies 
well within the phonon energies of the MgO crysta h 30 ' 31 
Therefore, rather than expecting to see a static Jalin- 
Teller distortion, as it is present in some Fe 2+ containing 
metal complexes,— ~— one should take this distortion as 
dynamical in which the phonons are strongly coupled to 
the electronic degrees of freedom as pointed out by Ham 
and co- workers! 22 ' 25 ' 35 To generate static distortions one 
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FIG. 6. (Color online) Simulations of the temperature dependence of the Fe L3 and L2 XAS spectra in the D^h symmetry. 
Parameters for the simulations are explained in the text. 







1 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 
O, Best fit D 

k 4h 








J=3 ^^^^ 










-""^75 -76 










L0- 


















... , 


. doublet 


5 ■ 


r * 














- 








■ 


> 1 












-20 -40 -60 -80 -100 

At 2g (meV) 



FIG. 7. (Color online) The total energy level diagram of the 
Fe 2+ Fe06 cluster as a function of the D^u crystal field split- 
ting expressed in terms of Af2 9 • 



would need an effective crystal field splitting of about 
0.2 eV or larger. 

It is now interesting to see what consequences the pres- 
ence of such a D^h crystal field splitting has for the mag- 
netic properties of the Fe 2+ ion. We calculate the spin 
and orbital contributions to the magnetic moments in 
the presence of an exchange field (H ex ) of 75 meV. We 
chose for this value as it may be taken as a crude es- 
timate for the case of FeaO,}. Fig. [5] shows the results 
of the calculations for several values of the crystal field 
energy At 2g . For At 2g — 0, i.e., the Oh case, the orbital 
moment is very large, very close to 1.0 (is- Yet, it also 
decreases rapidly with temperature: at 500 K it becomes 



0.5 /ib- This is the consequence of the thermal popula- 
tion of the J = 2 and 3 excited states. Upon switching 
on the D^h crystal field to -76 meV, the orbital moment 
gets also reduced, to about 0.7 (Lb already at K. In- 
creasing further the crystal field to -150 meV, -225 meV 
and -300 meV produces smaller and smaller orbital mo- 
ments, i.e., about 0.4, 0.3, and 0.2 (is, respectively. The 
spin moment, nevertheless, always stays close to about 

4 (JL B - 

These findings could provide an interesting path to 
critically test recent electronic structure theories^ for 
the explanation of the experimentally observed complex 
charge and orbital order phenomena in FeaO^ - — An ac- 
curate measurement of the orbital moment, i.e., the or- 
bital moment at the Fe 2+ sites (the Fe 3+ 3d 5 with their 
high-spin half-filled shell do not carry an orbital mo- 
ment), will provide a direct estimate of the magnitude 
of the effective crystal field splitting. This in turn will 
determine whether the occupied minority t 2g orbital of 
the Fe 2+ is made of mainly real space orbitals or has a 
more complex nature. Only for crystal fields substan- 
tially larger than the spin-orbit interaction one can ob- 
tain the real space orbitals necessary to build a robust or- 
bital ordering. In this sense the measurement of 0.76 (is 
orbital moment at 88 K by Huang et al~ would suggest 
the occupation of a complex t 2g orbital and a crystal field 
too small to produce a static Jahn- Teller distortion. On 
the other hand, the measurement by Goering et al£L of 
0.01 (Ib would support the scenario for real space orbitals 
and large static Jahn- Teller distortions, much more in 
line with the recent theoretical studies.— ^ Nevertheless, 
the issue on the magnitude of the orbital moment is not 
clear and is subject of debatei 38 ' 39 

To conclude, we have succeeded in preparing the 
Fe 2+ :MgO impurity system using MBE thin film technol- 
ogy. The resulting Fe L 2 ^ soft x-ray absorption spectra 
display very sharp features, thereby allowing us to firmly 
establish that the Fe local coordination has a lower sym- 
metry than Oh- Detailed analysis of the spectral line 
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FIG. 8. (Color online) Simulated temperature dependence of the spin (m s = 2-S x , panel (a)), and orbital (mi = L x , panel (b)) 
contributions to the local magnetic moment of the FeOe cluster for various D^h crystal field energies (At2 g ) and an exchange 
field of 75 meV. 



shape and its temperature dependence reveals that the 
local symmetry is D 4 h with an effective t 2g crystal field 
splitting of about -76 meV. With an energy well within 
the phonon frequencies of MgO, this gives rise to a dy- 
namic Jahn- Teller distortion. Using this Fe 2+ impurity 
system as a model we showed that an accurate measure- 
ment of the orbital moment in Fe304 will provide a di- 
rect estimate for the effective local low-symmetry crys- 
tal fields on the Fe 2+ sites, important for the theoretical 



modeling of the formation of orbital ordering. 
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